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I . INTRODUCTION 



A. NEED FOR STUDY 

Electromagnetic pulse has become a possible threat to 
nearly all sophisticated military systems [Ref. 1]. Under 
the proper circumstances a significant portion of the energy 
released during a nuclear detonation will be appear as an 
Electromagnetic Pulse (EMP) having the same frequencies as 
those employed by most of our commercial radio and military 
system equipments. 

Any metallic object, such as an aircraft structure, ex- 
posed to electromagnetic fields can be a collector of elec- 
tromagnetic energy, that is, act like an antenna, even 
though it was never intended to be that. 

The crossed-dipole receiving antenna has been used as 
a representative model to approximate electromagnetic pulse 
effects on aircraft [Ref. 2]. The crossed-dipole, in this 
case the aircraft, can be modeled in its simplest way as a 
two-dimensional structure , that is , a cross with thickness 
equal to zero . 

B. BACKGROUND INFORMATION 

Electromagnetic boundary-value problems can be solved 
exactly via classical separation of variable analysis only 
for a few geometries, such as the sphere, spheroid, circular 
cylinder, elliptical cylinder, strip and wedge. Such solu- 
tions can be expressed in terms of a summation of a set of 
eigenfunctions which can be evaluated with a high speed 
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computer. However, for problems involving complicated geo- 
metries the exact solution is not available. Thus, the ap- 
proximate and numerical methods are of great significance. 
The research discussed herein is directed toward the numeri 
cal solution of the problem of scattering from conducting 
bodies of arbitrary shape. 

Two methods [Ref. 3] are available for electromagnetic 
modeling of continuous conducting surfaces with arbitrary 
shape: the wire grid model and the surface current model 

using rectangular-pulse bases. Both methods have similar 
limitations with the maximum cell width restricted to ap- 
proximately X/10. Unless the conducting body is symmetric 
or is a figure of revolution, computer storage requirements 
have limited the moment-method application to bodies with 
surface area not exceeding one or two square wavelengths. 

In this paper, the piecewise-sinusoidal reaction tech- 
nique applied by Richmond [Ref. 3] to scattering by conduct 
ing bodies of arbitrary shape is used. 

C. STATEMENT OF THE PROBLEM 

The objective of this research is to find the current 
distribution in a two-dimensional cross. In order to get 
an understanding of the "junction problem" that will be 
present in the cross structure, the induced current density 
on a perfectly-conducting rectangular plate will be first 
solved for, this result having been presented by Richmond 

f 

and Wang [Ref. 3]. Likewise the square plate solved by 
Rahmat-Samii and Mitra [Ref. 4] will be studied. From this 
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the complexity of the structure will be increased, solving 
for current distribution on structures like an L, T and 
finally the cross . 

D. SCOPE AND LIMITATIONS OF THE STUDY 

The general solutions steps are : 

1) Divide the continuous conducting surface into cells 
with the surface current distribution being expanded in 
overlapping sinusoidal bases functions. 

2) Via an application of Galerkin’s method, the inte- 
gral equation formulated with the zero-reaction concept is 
reduced to a matrix equation. 

3) The current distribution over the conducting surface 
is determined via matrix inversion. 

The manner in which the structure rs subsectxoned is 
very important in any numerical analysis approach. Likewise 
the basic current distribution chosen for the subsections 
drastically affects the amount and difficulty of calculations. 
These methods will yield approximate results, the quality of 
the approximation depending on the amount of effort (and 
computer time) one is willing to devote to the problem. 

As will be seen symmetry may simplify considerably the 
number of calculations involved. In some cases one may want 
to take advantage of symmetry; in other cases it may be 
simpler to ignore it. It appears that one must trade some 
formulation time for computer time. Is the tradeoff worth- 
while? The answer to this question will depend on the extent 
of the computations. In reducing the number of unknowns, one 
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increases the analysis time but reduces computer programming 
time and computer running time. 
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II. METHOD OF ANALYSIS 



A. METHOD OF MOMENTS 

The method of moments has been widely explained in the 
literature. The interested reader may refer on this subject 
to References 5 , 6 and 7 . 

B. PIECEWISE SINUSOIDAL BASIS SETS 

The physical problem may be represented by many different 
operator equations, and a suitable one must be chosen. There 
are also an infinite number of sets of expansion functions 
and testing functions that may be chosen. 

Intuitively, the more nearly the expansion set resembles 
the actual form of the current distribution, the fewer the 
terms required to accurately represent the current. In ad- 
dition, certain types of basis functions may reduce many 
time-consuming numerical operations incorporated in approxi- 
mations. The piecewise-sinusoidal basis set exhibits these 
qualities. The objective is to represent the current dis- 
tribution, over individual subsections of a domain with ex- 
pansion functions that possess greater curve-fitting potential 
than do flat pulses and yet lead to well-conditioned matrices 
in approximate solution procedures. The application of piece- 
wise-sinusoids in Galerkin's method and the method of moments 
to electromagnetic field problems is known as the Piecewise- 
Sinusoidal Reaction Matching Technique. 
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C. PIECEWISE -SINUSOIDAL REACTION MATCHING TECHNIQUE. SURFACE MODEL 



The reaction concept and its application has been dis- 
cussed by Rumsey [Ref. 8], Cohen [Ref. 9], Harrington [Ref. 
10] and Richmond [Ref. 11]. 

The following material has been previously presented 
by Richmond and Wang [Ref. 3]; it is presented here for 
completeness . 

Consider the exterior scattering problem illustrated 
in Figure 1. In the presence of a dielectric or conducting 
body, the impressed electric and magnetic currents (J.pM.) 
generate the electric and magnetic field intensities (E ,H) . 
For simplicity let the exterior medium by free space. 

From the surface-equivalence theorem of Schelkunoff, 
the interior field will vanish (without disturbing the ex- 
terior field) if the following surface-current densities 
are introduced: 



J = n x H 



( 1 ) 



M = E x n 
— s — 



( 2 ) 



on the closed surface s of the scatterer. (The unit vector 

A 

n is directed outward on s . ) In this situation, illustrated 
in Figure 2, it is possible to replace the scatterer with 
free space without disturbing the field anywhere. 

By definition, the incident field (E^,H^) is generated 
by (J^>M^) in free space, and the scattered field is 



E 

— s 
H 



— s 




H - 



H. 



— l 



(3) 

(4) 
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Figure 1. 



The Source ( J . ,M. ) Generates the Field (E,H) with 
Scatterer . 




Figure 2. The Interior Field Vanishes V/ hen the Currents 
(J ,M ) are Introduced on the Surface of the 
Scatterer. 
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When the surface current (J ,M ) radiates in free space, it 
generates the field (E ,H ) in the exterior and (-E.,-H.) 
in the interior region. This result, illustrated in Figure 
3, is deduced from Figure 2 and the superposition theorem. 

With the scatterer replaced by free space, it can be 
noted in Figure 2 that the interior region has a null field. 
As shown in Figure 4, by placing an electric test source J 
in this region it is found from the reciprocity theorem 
that 



where (E^H^) is the free space field of the test source. 

In words, Equation (5) states that the interior test source 
has zero reaction with the other sources. This "zero reac- 
tion theorem" was developed by Rumsey [Ref. 8]. 

Equation (5) is the integral equation for the scattering 
problem, and the objective is to use this equation to determine 

the surface current distribution J and M . To accomplish 

s s 

this, these functions are expanded in finite series yielding 
a finite number N of unknown expansion constants . Next N 
simultaneous linear equations are obtained to permit a solu- 
tion for these constants. One such equation is obtained from 
Equation (5) each time a new test source is set up. 

If the scatterer is a perfect conductor, the magnetic 
current vanishes . Assuming a finite conductivity and 
using the impedance boundary condition 




(5) 




( 6 ) 
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Figure 3. The Exterior Scattered Field may be Generated by 

(J ,M ) in Free Space. 
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Figure 4. An Electric Test Source is Positioned in the 
Interior of the Scattering Region. 
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where Z_ denotes the surface impedance equation (5) yields 
- ff s J s *[E t -(nxH t )Z s ]ds = ///( J i *E t -M i *H t )dv. (7) 

For two-dimensional problems involving arbitrary scat- 

terers , J and M are functions only of the position on the 
— s —s 

surface of the scatterer. If vanishes and denoting 

(E ,H ) as the free-space field of test-source m, Equation 
m m 

(7) yields 

-//_ J • CE -(nxH )Z ]ds = ff J. -E ds . (8) 

s — s — m — m s l — m 



Now representing the electric current distribution as 
follows : 



J 



— s 





(9) 



where the complex constants I are samples of the function 

J . The vector functions J are known as basis functions. 

— s — n 

Basis functions J , and testing (weighting) sources with 

unit current density at their terminals are used. 

From Equations (8) and (9) the following simultaneous 
linear equation is obtained 



N 

I, I C = A 
n=l n mm m 



( 10 ) 



with m = 1,2,3,...,N, where 

C m =-// J .[E -(nxH )Z]ds 
mm — n — m s 

n 



-// J • E ds 
m —ra —n 



m 



ff. J. *E ds 
l — l -m 



ff J -E. ds 
m — m —l 



( 11 ) 

( 12 ) 
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The integrations in Equations (11) and (12) extend over 

the region where the integrand is non-zero. Region "m" 

covers the interior test source , J , and region "n" is that 

— m 

portion of the surface s covered by the expansion function 

J . The first and the second integrals in Equation (11) are 

related by the reciprocity theorem, where E^ is the free- 

space field generated by J n and the associated magnetic 

current M . 

— n 



D. UNIDIRECTIONAL CURRENT MODEL 

A two-dimensional surface dipole located on the y-z plane 
is illustrated in Figure 5(a). This source is an electric 
surface-current density with height 2D and width w. The 
surface current density is given by 



T _ sinK(D-z) 2 
-(z) " sinKD z 

T _ sinK(D+z) 

— (z) ~ sinKD 



for 0 < z < D (13) 

for-D < z < 0 (14) 



where K = 2 tt/X. 

As illustrated in Figure 5(b), the current density is as- 
sumed to be piecewise-sinusoidal in the z-direction, with 
period equal to X, vanishing at the edges z = D and z = -D, 
and equal to unity at the center. From Figure 5(c) it can 
be seen that the current density is assumed to be uniformly 
distributed in the tranverse direction. 

The dipole current distribution, Equations (13) and 

(14), will be used as the basis functions (J in Equation 

— n 

(9)) for expanding the unknown current distribution induced 
on a conducting surface. Also, surface dipoles will be 
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Figure 5. An Electric Surface Dipole and Its Current 
Density Distribution. 
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employed as test sources with the reaction concept to solve 
the integral equation. 

1 . Calculation of the Coupling Term 

The complex number C^ in Equation (11) represents 
the mutual coupling between two sources. Let <m,n> denote 
the mutual coupling between sources, then Equation (11) can 
be written as 

C = -<m,n> = -ff J *E ds (15) 

mn ’ m -mi n 



Relating the mutual coupling between two sources to the cir- 
cuit parameters [Ref. 3] 



C = V I 
mn mn on 



(16) 



where V is the open circuit terminal voltage induced at 
mn ° 

source m by source n, and I is the short circuit terminal 
current induced at source n. 

The surface dipole self impedance has been defined [Ref. 
3] with the induced-emf formulation: 



Z 



mm 



V 

mm 

I 

om 



I 



C 

mm 



I 



om om 



(17) 



From Equation (15), Equation (17) yields 



- 1 



mm 



I I 
om om 



ff J • E ds. 
m -mi -m 



(18) 



where J is the surface-current density of source m and E_ 
-m J — m 

is its free space electric field. 

The mutual impedance between two surface dipoles is 
defined [Ref. 3] by 
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